Iterative learning control method for multi-particle vehicle platoon driving system

ABSTRACT

The present invention discloses an iterative learning control (ILC) method for a multi-particle vehicle platoon driving system, and relates to the field of ILC. The method includes: firstly, discretizing a multi-particle train dynamic equation using a finite difference method to obtain a partial recurrence equation, and then transforming the partial recurrence equation into a spatially interconnected system model; secondly, transforming the spatially interconnected system model into an equivalent one-dimensional dynamic model using a lifting technology, and in order to compensate input delay, designing an ILC law based on a state observer, and thirdly, transforming a controlled object into an equivalent discrete repetitive process according to the ILC law, and converting a controller combination problem into a linear matrix inequality based on stability analysis of the repetitive process. The method is simple and easy to implement, considers structure uncertainty of the system, and has a good control performance and robustness.

CROSS-REFERENCE TO RELATED APPLICATIONS

The application claims priority to Chinese patent application No. 202110680017.6, filed on Jun. 18, 2021, the entire contents of which are incorporated herein by reference.

TECHNICAL FIELD

The present invention relates to the field of iterative learning control (ILC), and particularly to an ILC method for a multi-particle vehicle platoon driving system, which is used for multi-particle vehicle platoon driving systems with structure uncertainty and input delay.

BACKGROUND

In recent years, with a rapid development of railways, an efficient travel mode of trains is popular with people. A multi-particle vehicle platoon driving system is a spatially interconnected system formed by coupling a plurality of train carriages to each other and has characteristics of time-space coupling, a large number of variables and a high dimension. On the premise of guaranteeing safe operation of the trains, an improvement on an operation speed and performance of the trains in limited time and trials becomes a focus. However, due to a complex operation environment of the trains, the system always has uncertainty in practice. Furthermore, information interaction transmission exists in control implementation, and input delay is usually involved, resulting in delay in output response of the system to affect stability and a performance of the system. In view of operation characteristics of the multi-particle vehicle platoon driving system, ILC is used to solve a tracking problem of the multi-particle vehicle platoon driving system.

ILC is a new intelligent control method which can improve a tracking performance of the controlled system in limited time by repeatedly learning from previous trials. ILC has most remarkable characteristics of a small calculated amount, low requirements for prior knowledge of system dynamics, simplicity and implementation easiness. In the present application, application of ILC to the multi-particle vehicle platoon driving system with uncertainty and input delay may increase a response speed and trajectory tracking precision of the system and improve a control performance.

SUMMARY

In view of the above problems and technical demands, the inventor proposes an ILC method for a multi-particle vehicle platoon driving system. The present invention has the following technical solution.

An ILC method for a multi-particle vehicle platoon driving system includes the following steps:

step 1: establishing a spatially interconnected system model of the multi-particle vehicle platoon driving system;

wherein a dynamic equation of the multi-particle vehicle platoon driving system is described as

m{umlaut over (x)} _(i)({tilde over (t)})=k(x _(i+1)({tilde over (t)})−x _(i)({tilde over (t)}))−k(x _(i)({tilde over (t)})−x _(i−1)({tilde over (t)}))+u _(i)({tilde over (t)}−τ) y _(i)({tilde over (t)})=x _(i)({tilde over (t)})  (1)

wherein x_(i)({tilde over (t)}) represents a position of an ith train, u_(i)({tilde over (t)}) is control input and represents traction of the train, τ is a time-delay constant caused by signal transmission, y_(i)({tilde over (t)}) is measurement output and represents a speed of the train, m denotes mass of each train and has a variation range±ω _(m), and k is a spring coefficient and has a variation range±ω _(k);

sampling time T is selected, and the equation (1) is approximately discretized using a finite difference method; that is,

${\left( {{\overset{˙}{x}}_{i}\left( \overset{\sim}{t} \right)} \right)_{t,s} = \frac{{x_{1}\left( {{t + 1},s} \right)} - {x_{1}\left( {t,s} \right)}}{T}},$ $\left( {{\overset{¨}{x}}_{i}\left( \overset{\sim}{t} \right)} \right)_{t,s} = \frac{{x_{2}\left( {{t + 1},s} \right)} - {x_{2}\left( {t,s} \right)}}{T}$

wherein t and s are discrete time and a train serial number respectively, and substitution of the above formula into the equation (1) gives a partial recurrence equation:

$\begin{matrix} {{x_{1}\left( {{t + 1},s} \right)} = {{x_{1}\left( {t,s} \right)} + {T{x_{2}\left( {t,s} \right)}}}} & (2) \end{matrix}$ ${x_{2}\left( {{t + 1},s} \right)} = {{{- 2}\frac{kT}{m}{x_{1}\left( {t,s} \right)}} + {x_{2}\left( {t,s} \right)} + {\frac{kT}{m}{x_{1}\left( {t,{s + 1}} \right)}} + {\frac{kT}{m}{x_{1}\left( {t,{s - 1}} \right)}} + {\frac{T}{m}{u\left( {{t - \tau},s} \right)}}}$ y(t, s) = x₂(t, s)

assuming that information transferred between the trains is their position information, that is, interconnected variables are set as: w+(t,s)=w.(t,s)=x₁(t,s), v+(t,s)=x₁(t,s−1), and v(t,s)⁻x₁(t,s+1), and an output variable is set as: y(t,s)=x₂(t,s), the equation (2) is transformed into an uncertain spatially interconnected system model, i.e.,

$\begin{matrix} {\begin{bmatrix} \begin{matrix} \begin{matrix} {x\left( {{t + 1},s} \right)} \\ {w\left( {t,s} \right)} \end{matrix} \\ {q\left( {t,s} \right)} \end{matrix} \\ {y\left( {t,s} \right)} \end{bmatrix} = {\begin{bmatrix} {A_{11}\text{?}} & {A_{12}\text{?}} & {B_{11}\text{?}} & {B_{11}\text{?}} \\ {A_{21}\text{?}} & 0 & 0 & 0 \\ {C_{11}\text{?}} & {C_{12}\text{?}} & {D_{11}\text{?}} & {D_{12}\text{?}} \\ {C_{21}\text{?}} & 0 & 0 & 0 \end{bmatrix}\begin{bmatrix} \begin{matrix} \begin{matrix} {x\left( {t,s} \right)} \\ {v\left( {t,s} \right)} \end{matrix} \\ {p\left( {t,s} \right)} \end{matrix} \\ {u\left( {{t - \tau},s} \right)} \end{bmatrix}}} & (3) \end{matrix}$ ?indicates text missing or illegible when filed

wherein

${{A_{11}\text{?}} = \begin{bmatrix} 1 & T \\ {{- 2}\frac{kT}{m}} & 1 \end{bmatrix}},{{A_{12}\text{?}} = \begin{bmatrix} 0 & 0 \\ \frac{kT}{m} & \frac{kT}{m} \end{bmatrix}},{A_{21}{\text{?}\begin{bmatrix} 1 & 0 \\ 1 & 0 \end{bmatrix}}},$ ${{B_{11}\text{?}} = \begin{bmatrix} 0 & 0 \\ \frac{T}{m} & {- \frac{T}{m}} \end{bmatrix}},{{B_{12}\text{?}} = \begin{bmatrix} 0 \\ \frac{T}{m} \end{bmatrix}},{{C_{11}\text{?}} = \begin{bmatrix} {{- 2}{\overset{\_}{\omega}}_{k}} & 0 \\ {{- 2}\frac{{\overset{\_}{\omega}}_{m}k}{m}} & 0 \end{bmatrix}},$ ${{C_{12}\text{?}} = \begin{bmatrix} {\overset{\_}{\omega}}_{k} & {\overset{\_}{\omega}}_{k} \\ \frac{{\overset{\_}{\omega}}_{m}k}{m} & \frac{{\overset{\_}{\omega}}_{m}k}{m} \end{bmatrix}},{{C_{21}\text{?}} = {{\begin{matrix} \left\lbrack 0 \right. & {\left. 1 \right\rbrack,} \end{matrix}D_{11}\text{?}} = \begin{bmatrix} 0 & 0 \\ \frac{{\overset{\_}{\omega}}_{m}}{m} & {- \frac{{\overset{\_}{\omega}}_{m}}{m}} \end{bmatrix}}},$ ${D_{12}\text{?}} = \begin{bmatrix} 0 \\ \frac{{\overset{\_}{\omega}}_{m}}{m} \end{bmatrix}$ ?indicates text missing or illegible when filed

in the above formula, x(t,s)∈ R^(n), u(t,s)∈ R^(q) and y(t,s)∈ R represent state, input and output variables of an sth subsystem respectively; τ(0<τ<α) denotes the time-delay constant; w(t,s)=[w₊ ^(T)(t,s)w^(T)(t,s)]^(T) and v(t,S)=[v₊ ^(T)(t,S) v⁻ ^(T)(t,S)]^(T) represent interconnection between adjacent subsystems, and satisfy

w ₊(s)=v ₊(s+1) w(s)=v(s−1)  (4)

boundary conditions thereof are set as: v₊(1)=w.(1)=0 and v.(n)=w₊(n)=0, and n is a number of the subsystems;

p(t,s) denotes a pseudo-input channel of structure uncertainty, q(t,s) denotes a pseudo-output channel of the structure uncertainty, and

p(t,s)θ_(s) q(t,s)  (5)

wherein an uncertainty block is: θs=diag{δ₁I_(r), . . . , δ₁I_(r)}, δ1 ∈0R, |δ_(i)|≤1. i=1, . . . , f; δ_(i) describes a variation of a dynamic parameter of the system, and I_(r) _(i) is a unit matrix with a dimension of r_(f);

step 2: transforming the spatially interconnected system model;

wherein the model (3) is transformed into an equivalent one-dimensional dynamic model using a lifting technology, and lifting vectors are defined as follows:

-   -   X(t)=[x(t,1)^(T), x(t,2)^(T), . . . , x(t,n)^(T)]^(T)     -   V(t)==[v(t,1)^(T), v(t,2)^(T), . . . , v(t,n)^(T)]^(T)     -   W(t)=[w(t,1)^(T), w(t,2)^(T), . . . , w(t,n)^(T)]^(T)     -   P(t)=[p(t,1)^(T), p(t,2)^(T), . . . , p(t,n)^(T)]^(T)     -   Q(t)=[q(t,1)^(T), q(t,2)^(T), . . . , q(t,n)^(T)]^(T)     -   U(t−τ)=[u(t−τ, 1)^(T),u(t−τ,2)^(T), . . . , u(t−τ,n)^(T)]^(T)     -   Y(t)=[y(t,1)^(T), y(t,2)^(T), . . . , y(t,n)^(T)]^(T)

then, the whole uncertain spatially interconnected system model can be equivalently described by the following model:

$\begin{matrix} {\begin{bmatrix} {X\left( {t + 1} \right)} \\ {W(t)} \\ {Q(t)} \\ {Y(t)} \end{bmatrix} = {\begin{bmatrix} A_{11} & A_{12} & B_{11} & B_{12} \\ A_{21} & 0 & 0 & 0 \\ C_{11} & C_{12} & D_{11} & D_{12} \\ C_{21} & 0 & 0 & 0 \end{bmatrix}\begin{bmatrix} {X(t)} \\ {V(t)} \\ {P(t)} \\ {U\left( {t - \tau} \right)} \end{bmatrix}}} & (6) \end{matrix}$ whereinA₁₁ = diag{A₁₁¹, …, A₁₁^(n)}, A₁₂ = diag{A₁₂¹, …, A₁₂^(n)}, A₂₁ = diag{A₂₁¹, …, A₂₁^(n)}, B₁₁ = diag{B₁₁¹, …, B₁₁^(n)}, B₁₂ = diag{B₁₂¹, …, B₁₂^(n)}, C₁₁ = diag{C₁₁¹, …, C₁₁^(n)} C₁₂ = diag{C₁₂¹, …, C₁₂^(n)}, C₂₁ = diag{C₂₁¹, …, C₂₁^(n)}, D₁₁ = diag{D₁₁¹, …, D₁₁^(n)}, D₁₂ = diag{D₁₂¹, …, D₁₂^(n)},

the model (6) contains the interconnected variables, and therefore, the model (6) is required to be further simplified;

based on an interconnection characteristic of the formula (4) and the boundary conditions thereof, an equation relationship between the interconnected variables is obtained:

W(t)=ηV(t)  (7)

wherein q is a permutation matrix unrelated to time t;

the formula (7) is substituted into (6) to obtain

V(t)=η⁻¹ A ₂₁ X(t)  (8)

and then, the formula (8) is substituted into (6) to eliminate the interconnected variables W(t) and V(t), so as to obtain the following equivalent uncertain model:

$\begin{matrix} {\begin{bmatrix} \begin{matrix} {X\left( {t + 1} \right)} \\ {Q(t)} \end{matrix} \\ {Y(t)} \end{bmatrix} = {\begin{bmatrix} {\overset{\_}{A}}_{11} & B_{11} & B_{12} \\ {\overset{\_}{C}}_{11} & D_{11} & D_{12} \\ C_{21} & 0 & 0 \end{bmatrix}\begin{bmatrix} \begin{matrix} {X(t)} \\ {P(t)} \end{matrix} \\ {U\left( {t - \tau} \right)} \end{bmatrix}}} & (9) \end{matrix}$

wherein

Ā ₁₁ =A ₁₁ +A ₁₂η¹ A ₂₁

C ₁₁ =C ₁₁ +C ₁₂η⁻¹ A ₂₁

according to the formula (5), the following formula is obtained:

P(t)=θQ(t)  (10)

wherein the uncertainty block is: θ=diag{θ₁, . . . ,θ_(n)}, θ_(i)≤I, i=1, . . . ,n;

the formula (10) is substituted into (9) to obtain

Q(t)=(1−D ₁₁θ)¹( C ₁₁ X(t)+D ₁₂ U(t−τ))  (11)

then, the formula (11) is substituted into (9) to eliminate uncertainty variables P(t) and Q(t) using an elimination method, so as to obtain a general-form state space model:

X(t+1)=(A+ΔA)X(t)+(B+ΔB)U(t−τ) Y(t)=CX(t)  (12)

wherein

-   -   A=Ā₁₁, B=B₁₂, C=C₂₁     -   ΔA=B₁₁θ(I−D₁₁θ)⁻¹ C ₁₁     -   ΔB=B₁₁θ(I−D₁₁θ)⁻¹D₁₂

step 3: designing an ILC law based on a state observer,

wherein the state space model (12) is described as an ILC structure form:

X _(k+1)(t+1)=(A+ΔA)X _(k+1)(t)+(B+ΔB)U _(k+1)(t−τ) Y _(k+1)(t)=CX _(k+1)(t)  (13)

wherein k+1 represents a current operation trial of the system, t ∈[0,α] denotes a limited working cycle of each trial of the system, and the input-delay constant satisfies a condition: τ<α;

then, the ILC law is expressed as

U _(k+1)(t)=U _(k)(t)+r _(k+1)(t)  (14)

the current control signal U_(k+1)(t) is equal to the control signal U_(k)(t) of a previous trial plus the update term r_(k+1)(t) which is calculated from previous error information;

a system tracking error in a (k+1)th trial is

e _(k−1)(t)=Y _(r)(t−τ)−Y _(k+1)(t)  (15)

wherein Y_(r)(t) is a desired output trajectory;

considering delay in output response, the tracking error is re-described as

e _(k−1)(t)=Y _(r)(t−τ)−Y _(k−1)(t)  (16)

a state error vector is introduced:

X _(k+1)(t+1)=X _(k+1)(t)−X _(k)(t)  (17)

assuming that Y_(r)(0)=Y_(k)(0)=CX_(k)(0), and X _(k) (0)=0, that is, the system is returned to a same initial state in each trial,

$\begin{matrix} \begin{matrix} {{{\overset{\_}{X}}_{k + 1}\left( {t + 1} \right)} = {{X_{k + 1}(t)} - {X_{k}(t)}}} \\ {= {{\left( {A + {\Delta A}} \right)\left\lbrack {{X_{k + 1}\left( {t - 1} \right)} - {X_{k}\left( {t - 1} \right)}} \right\rbrack} +}} \\ {\left( {B + {\Delta B}} \right)\left\lbrack {{U_{k + 1}\left( {t - 1 - \tau} \right)} - {U_{k}\left( {t - 1 - \tau} \right)}} \right\rbrack} \\ {= {{\left( {A + {\Delta A}} \right){{\overset{\_}{X}}_{k + 1}(t)}} + {\left( {B + {\Delta B}} \right){r_{k + 1}\left( {t - 1 - \tau} \right)}}}} \end{matrix} & (18) \end{matrix}$ and $\begin{matrix} \begin{matrix} {{e_{k + 1}(t)} = {{Y_{r}\left( {t - \tau} \right)} - {Y_{k + 1}(t)}}} \\ {= {{e_{k}(t)} - {C\left\lbrack {{X_{k + 1}(t)} - {X_{k}(t)}} \right\rbrack}}} \\ {= {{{- {C\left( {A + {\Delta A}} \right)}}{{\overset{\_}{X}}_{k + 1}(t)}} - {{C\left( {B + {\Delta B}} \right)}{r_{k + 1}\left( {t - 1 - \tau} \right)}} + {e_{k}(t)}}} \end{matrix} & (19) \end{matrix}$

in order to compensate input delay, the following state observer is constructed based on output information in the current trial:

{tilde over (X)} _(k+1)(t+1)=A{tilde over (X)} _(k+1)(t)+Br _(k+1)(t−1)+L[ Y _(k+1)(t)− Y _(k+1)(t−τ)]  (20)

wherein Y _(k+1)(t)=Y_(k+1)(t−1)−Y_(k)(t−1)=CX _(k+1)(t) and {tilde over (X)}_(k+1)(t) is a t-step ahead prediction of the state X _(k+1)(t); that is, X_(k+1)(t) is an estimation value of X _(k+1)(t+τ); and L is observer gain to be designed;

an observation error is defined as

{tilde over (e)} _(k+1)(t)= X _(k+1)(t)−{tilde over (X)} _(k+1)(t−τ)  (21)

it is assumed that the update term in the ILC law (14) is

r _(k+1)(r)=K ₁ {tilde over (X)} _(k+1)(t+1)+K ₂ e _(k)(t+τ)+K ₃ e _(k)(t+1+τ)−e _(k)(t+τ))  (22)

wherein K₁, K₂ and K₃ are learning gain to be designed;

the update term is formed by state feedback information and PD-type previous tracking error information, and when the learning gain K₂ is equal to K₃, the formula (22) is simplified to P-type ILC;

(22) is substituted into (20) to obtain

$\begin{matrix} \begin{matrix} {{{\overset{\sim}{X}}_{k + 1}\left( {t + 1 - \tau} \right)} = {{\left( {A + {BK}_{1}} \right){{\overset{\sim}{X}}_{k + 1}\left( {t - \tau} \right)}} + {{LC}{{\overset{\sim}{e}}_{k + 1}\left( {t - \tau} \right)}} +}} \\ {{{B\left( {K_{2} - K_{3}} \right)}{e_{k}\left( {t - 1} \right)}} + {{BK}\text{?}{e_{k}(t)}}} \end{matrix} & (23) \end{matrix}$ and $\begin{matrix} \begin{matrix} {{{\overset{\sim}{e}}_{k + 1}\left( {t + 1} \right)} = {{{\overset{\_}{X}}_{k + 1}\left( {t + 1} \right)} - {{\overset{\sim}{X}}_{k + 1}\left( {t + 1 - \tau} \right)}}} \\ {= {{\left( {{\Delta A} + {\Delta{BK}_{1}}} \right){{\overset{\sim}{X}}_{k + 1}\left( {t - \tau} \right)}} + {\left( {A + {\Delta A}} \right){{\overset{\sim}{e}}_{k + 1}(t)}} -}} \\ {{{LC}{{\overset{\sim}{e}}_{k + 1}\left( {t - \tau} \right)}} + {\Delta{B\left( {K_{2} - K_{3}} \right)}{e_{k}\left( {t - 1} \right)}} + {\Delta{BK}_{3}{e_{k}(t)}}} \end{matrix} & (24) \end{matrix}$ ?indicates text missing or illegible when filed

it is assumed that ζ_(k+1)(t)=[{tilde over (X)}_(k+1) ¹(t−τ){tilde over (e)}_(k+1) ⁻¹(t)]¹, and K=K₂-K₃, so as to obtain the following linear discrete repetitive process model with time delay:

$\begin{matrix} {{\zeta_{k + 1}\left( {t + 1} \right)} = {{\hat{A}\zeta\text{?}(t)} + {\hat{A}\text{?}{\zeta_{k + 1}\left( {t - \tau} \right)}} + {{\hat{B}}_{1}e\text{?}\left( {t - 1} \right)} + {{\hat{B}}_{0}e\text{?}(t)}}} & (25) \end{matrix}$ e_(k + 1)(t) = Ĉζ?(t) + D̂₁e_(k)(t − 1) + D̂₀e_(k)(t) wherein ${\hat{A} = \begin{bmatrix} {A + {BK}_{1}} & 0 \\ {{\Delta A} + {\Delta{BK}_{1}}} & {A + {\Delta A}} \end{bmatrix}},{{\hat{A}}_{t} = \begin{bmatrix} 0 & {LC} \\ 0 & {- {LC}} \end{bmatrix}}$ ${{\hat{B}}_{1} = \begin{bmatrix} {BK} \\ {\Delta{BK}} \end{bmatrix}},{{\hat{B}}_{0} = \begin{bmatrix} {BK}_{3} \\ {\Delta{BK}_{3}} \end{bmatrix}}$ Ĉ = [−C(A + ΔA) − C(B + ΔB)K₁ − C(A + ΔA)] D̂₁ = −C(B + ΔB)K, D̂₀ = I − C(B + ΔB)K₃ ?indicates text missing or illegible when filed

step 4: performing system stability analysis and learning gain solving on the linear discrete repetitive process model;

wherein a Lyapunov function is selected as

$\begin{matrix} {{V\left( {k,t} \right)} = {{V_{1}\left( {t,k} \right)} + {V_{2}\left( {k,t} \right)}}} & (26) \end{matrix}$ ${V_{1}\left( {t,k} \right)} = {{{\zeta_{k + 1}^{T}(t)}S{\zeta_{k + 1}(t)}} + {\sum\limits_{i = 1}{\text{?}{\zeta_{k + 1}^{T}\left( {t - i} \right)}Q{\zeta_{k + 1}\left( {t - i} \right)}}}}$ V₂(k, t) = e_(k)^(T)(t − 1)P₁e_(k)(t − 1) + e^(T)?(t)(P₂ − P₁)e?(t) whereinS = diag{S₁, S₂} > 0, Q = diag{Q₁, Q₂} > 0, andP₂ > P₁ > 0; ?indicates text missing or illegible when filed

V₁ (t,k) represents an energy change along a trial, and V₂(k,t) represents an energy change between the trials; and increments of subfunctions are

$\begin{matrix} {{\Delta{V_{1}\left( {t,k} \right)}} = {{\zeta_{k + 1}\text{?}\left( {t + 1} \right)S{\zeta_{k + 1}\left( {t + 1} \right)}} - {{\zeta_{k + 1}^{T}(t)}\left( {S - Q} \right){\zeta_{k + 1}(t)}} - {{\zeta_{k + 1}^{T}\left( {t - \tau} \right)}Q\zeta\text{?}\left( {t - \tau} \right)}}} & (27) \end{matrix}$  = H_(k + 1)^(T)(t)Φ₁^(T)SΦ₁H_(k + 1)(t) − ζ_(k + 1)^(T)(t)(S − Q)ζ_(k + 1)(t) − ζ_(k + 1)^(T)(t − τ)Qζ_(k + 1)(t − τ) $\begin{matrix} {{\Delta{V_{2}\left( {k,t} \right)}} = {{{e_{k + 1}^{T}(t)}P_{2}{e_{k + 1}(t)}} - {e^{T}\text{?}\left( {t - 1} \right)P_{1}e\text{?}\left( {t - 1} \right)} - {e^{T}\text{?}(t)\left( {P_{2} - P_{1}} \right)e\text{?}(t)}}} & (28) \end{matrix}$  = H_(k + 1)^(T)(t)Φ₂?P₂Φ₂H_(k + 1)(t) − e_(k)?(t − 1)P₁e?(t − 1) − e_(k)^(T)(t)(P₂ − P₁)e_(k)(t) wherein ${{H_{k + 1}(t)} = \begin{bmatrix} {\zeta_{k + 1}(t)} \\ {\zeta_{k + 1}\left( {t - \tau} \right)} \\ {e_{k}\left( {t - 1} \right)} \\ {e_{k}(t)} \end{bmatrix}},{\Phi_{1} = \begin{bmatrix} \hat{A} & {\hat{A}}_{\tau} & {\hat{B}}_{1} & {\hat{B}}_{0} \end{bmatrix}},{\Phi_{2} = \begin{bmatrix} \hat{C} & 0 & {\hat{D}}_{1} & {\hat{D}}_{0} \end{bmatrix}}$ ?indicates text missing or illegible when filed

a total increment of the function is

$\begin{matrix} {{\Delta{V\left( {k,t} \right)}} = {{{\Delta{V_{1}\left( {t,k} \right)}} + {\Delta{V_{2}\left( {k,t} \right)}}} = {{H_{k + 1}^{T}(t)}{\prod{H_{k + 1}(t)}}}}} & (29) \end{matrix}$ wherein $\prod{= {{{\begin{bmatrix} \hat{A} & {\hat{A}}_{\tau} & {\hat{B}}_{1} & {\hat{B}}_{0} \\ \hat{C} & 0 & {\hat{D}}_{1} & {\hat{D}}_{0} \end{bmatrix}^{T}\begin{bmatrix} S & 0 \\ * & P_{2} \end{bmatrix}}\begin{bmatrix} \hat{A} & {\hat{A}}_{\tau} & {\hat{B}}_{1} & {\hat{B}}_{0} \\ \hat{C} & 0 & {\hat{D}}_{1} & {\hat{D}}_{0} \end{bmatrix}} - \begin{bmatrix} {S - Q} & 0 & 0 & 0 \\ * & Q & 0 & 0 \\ * & * & P_{1} & 0 \\ * & * & * & {P_{2} - P_{2}} \end{bmatrix}}}$

if ΔV(k,t)<0 is satisfied for all k and t, the model (25) is stable along the trials, which has an equivalent condition of Π<0; a Schur complement lemma is applied to Π<0, left and right sides of an inequality are multiplied by diag{I, I, S¹, S¹, P₂ ¹, P₂ ¹}, variable substitution is performed, and it is assumed that S₁ ¹=W₁, S₂ ¹=W₂, S₁ ¹Q₁S₁ ¹=X₁, S₂ ¹Q₂S₂ ¹=X₂, P₂ ¹P₁P₂ ¹=Z₁, P₂ ¹=Z₂, K₁S₁ ¹=K₁W₁=R₁, LCS₂ ¹=LCW₂=R, KP₂ ¹=KZ₂=R₂, K₂P₂ ¹=K₁Z₂R₃, so as to obtain the following conclusion:

for the nominal linear discrete repetitive process model with time delay described in the formula (25), if there exist matrices W=diag(W₁,W₂)>0, X=diag{X₁,X₂}>0, Z₁>0, Z₂>0 and R, R₁, R₂ and R₃ satisfying the following linear matrix inequality:

$\begin{matrix} {\Xi_{1} = \text{ }{\begin{bmatrix} {- W_{1}} & 0 & 0 & {{AW}_{1} + {BR}_{1}} & 0 & 0 & R & {BR}_{2} & {BR}_{1} \\ * & {- W_{2}} & 0 & 0 & {AW}_{2} & 0 & {- R} & 0 & 0 \\ * & * & {- Z\text{?}} & {- \text{?}} & {- {CAW}_{2}} & 0 & 0 & {- {CBR}\text{?}} & {{Z\text{?}} - {{CBR}\text{?}}} \\ * & * & * & {{- W_{1}} + X_{1}} & 0 & 0 & 0 & 0 & 0 \\ * & * & * & * & {{- W_{2}} + X_{2}} & 0 & 0 & 0 & 0 \\ * & * & * & * & * & {- X_{1}} & 0 & 0 & 0 \\ * & * & * & * & * & * & X_{2} & 0 & 0 \\ * & * & * & * & * & * & * & {- Z_{1}} & 0 \\ * & * & * & * & * & * & * & * & {{- Z_{2}} + Z_{1}} \end{bmatrix} < 0}} & (30) \end{matrix}$ ?indicates text missing or illegible when filed

the model (25) is stable along the trials, and the learning gain of the update term (22) and the gain of the state observer (20) are:

K ₁ −R ₁ W ₁ ¹ ,K ₂ =R ₂ Z ₂ ¹ +R ₁ Z ₂ ¹ ,K ₃ =R ₁ Z ₂ ¹ ,L=RW ₂ ¹ C ¹(CC ¹)¹  (31)

when the structure uncertainty of the system is considered, the Schur complement lemma is applied to Π<0, the left and right sides of the inequality are multiplied by diag|I,I, S₁, S₁, P₂ ¹, P₂ ¹|, variable substitution is performed, and it is assumed that S₁ ¹=W₁, S₂ ¹=W₂, S₁ ¹Q₁S₁ ¹=X₁, S₂ ¹Q₂S₂ ¹=X₂, P₂ ¹P₁P₂ ¹=Z₁, P₂ ¹=Z₂, K₁S₁ ¹=K₁W₁R₁, LCS₂ ¹=LCW₂=R, KP₂ ¹=KZ₂=R₂, K₃P₂ ¹=K₃Z₂=R₃, so as to obtain

Ξ₁ +MΘN+N ^(T)Θ^(T) M ^(T)<0  (32)

wherein

-   -   M=[0 B₁₁ ^(T)−B₁₁ ^(T)C^(T) 0 0 0 0 0 0]^(r)     -   N=[0 0 0 C ₁₁W₁+D₁₂R₁ C ₁₁W₂ 0 0 D₁₂R₂ D₁₂R₃]     -   Θ=θ(I−D₁₁←)⁻¹, θ^(T)θ≤I

the formula (32) is equivalent to

$\begin{matrix} {{\Xi_{1} + {{\begin{bmatrix} \varepsilon^{- 1} & N^{T} & {\varepsilon M} \end{bmatrix}\begin{bmatrix} I & {- D_{11}} \\ * & I \end{bmatrix}}^{- 1}\begin{bmatrix} {\varepsilon^{- 1}N} \\ {\varepsilon M^{T}} \end{bmatrix}}} < 0} & (33) \end{matrix}$

wherein ε>0;

according to the Schur complement lemma, the formula (33) is described as

$\begin{bmatrix} \Xi_{1} & {\varepsilon^{- 1}N^{T}} & {\varepsilon M} \\ * & {- I} & D_{11} \\ * & * & {- I} \end{bmatrix} < 0$

left and right sides of the above formula are multiplied by diag(I,εI, εI), and ε² is replaced by ε, so as to obtain the following conclusion:

for the uncertain linear discrete repetitive process model with time delay described in the formula (25), if there exist the matrices W=diag{W₁,W₂}>0, X=diag{X₁,X₂}>0. Z₁>0, Z₂>0 and R, R₁, R₂ and R₃ satisfying the following linear matrix inequality:

$\begin{matrix} {\begin{bmatrix} \Xi_{1} & \Xi_{2} \\ * & \Xi_{3} \end{bmatrix} < 0} & (34) \end{matrix}$ wherein $\Xi_{2} = \begin{bmatrix} 0 & 0 & 0 & {{{\overset{\_}{C}}_{11}W_{1}} + {D_{12}R_{1}}} & {{\overset{\_}{C}}_{11}W_{2}} & 0 & 0 & {D_{12}R_{2}} & {D_{12}R_{3}} \\ 0 & {\varepsilon B_{11}^{T}} & {- \varepsilon B_{11}^{T}C^{T}} & 0 & 0 & 0 & 0 & 0 & 0 \end{bmatrix}^{T}$ $\Xi_{3} = \begin{bmatrix} {- \varepsilon I} & {\varepsilon D_{11}} \\ * & {- \varepsilon I} \end{bmatrix}$

the model (25) is robustly stable along the trials, and the learning gain of the update term (22) and the gain of the state observer (20) are given by the formula (31).

The present invention has the following beneficial technical effects.

In the present application, the robust ILC method for a multi-particle vehicle platoon driving system with structure uncertainty and input delay is researched; firstly, the multi-particle train dynamic equation is discretized using the finite difference method to obtain the partial recurrence equation, and then, the partial recurrence equation is transformed into the spatially interconnected system model; next, the spatially interconnected system model is transformed into the equivalent one-dimensional dynamic model using the lifting technology, and in order to compensate the input delay, the ILC law is designed based on the state observer, a controlled object is transformed into an equivalent discrete repetitive process according to the designed ILC law, sufficient conditions for the stability of the system along the trials are obtained according to a repetitive process stability theory, and a controller combination problem is converted into the linear matrix inequality. The method is simple and easy to implement, solves the problem of the delay in the output response caused by the input delay, considers the structure uncertainty of the system, and has a good control performance and robustness.

BRIEF DESCRIPTION OF DRAWINGS

FIG. 1 is a structural diagram of a multi-particle vehicle platoon driving system in the present application;

FIG. 2 is a structural diagram of a spatially interconnected system model in the present application;

FIG. 3 is an output curve of a first train particle under a nominal situation in the present application;

FIG. 4 is an output curve of a second train particle under the nominal situation in the present application;

FIG. 5 is an output curve of a third train particle under the nominal situation in the present application;

FIG. 6 is an RMS (Root-Mean-Square) contrast curve of a nominal system in the present application;

FIG. 7 is an output curve of the first train particle under an uncertain situation in the present application;

FIG. 8 is an output curve of the second train particle under the uncertain situation in the present application;

FIG. 9 is an output curve of the third train particle under the uncertain situation in the present application; and

FIG. 10 is an RMS contrast curve of an uncertain system in the present application.

DESCRIPTION OF EMBODIMENTS

An embodiment of the present invention is further described below with reference to the drawings.

FIG. 1 is a structural diagram of a multi-particle vehicle platoon driving system.

For the dynamic equation (1) of the system, the time-delay constant is τ=4; the mass of each train is m=2 [kg], and the variation range is ω _(m)=0.2 [kg]; the spring coefficient is k=2 [N/m], and the variation range is ω _(k)=0.2 [N/m].a the sampling time is T=0.05[s].

FIG. 2 is a structural diagram of a spatially interconnected system model, and all parameter matrices of the model (3) are

${A_{11}^{s} = \begin{bmatrix} 1 & 0.05 \\ {- 0.1} & 1 \end{bmatrix}},{A_{12}^{s} = \begin{bmatrix} 0 & 0 \\ 0.05 & 0.05 \end{bmatrix}},{A_{21}^{s} = \begin{bmatrix} 1 & 0 \\ 1 & 0 \end{bmatrix}},{B_{11}^{s} = \begin{bmatrix} 0 & 0 \\ 0.025 & 0.025 \end{bmatrix}}$ ${{B_{12}\text{?}} = \begin{bmatrix} 0 \\ 0.025 \end{bmatrix}},{{C_{11}\text{?}} = \begin{bmatrix} {- 0.4} & 0 \\ {- 0.4} & 0 \end{bmatrix}},{{C_{12}\text{?}} = \begin{bmatrix} 0.2 & 0.2 \\ 0.2 & 0.2 \end{bmatrix}},{{C_{21}\text{?}} = \begin{bmatrix} 0 & 1 \end{bmatrix}}$ ${D_{11}^{s} = \begin{bmatrix} 0 & 0 \\ 0.1 & {- 0.1} \end{bmatrix}},{D_{12}^{s} = \begin{bmatrix} 0 \\ 0.1 \end{bmatrix}}$ ?indicates text missing or illegible when filed

Considering a multi-particle vehicle platoon driving system with three interconnected trains, and supposing that an initial condition of the system is x₀(0,1)=x₀(0,2)=x₀(0,3)=0, and a limited working cycle in each trial is 20s, a speed reference trajectory of each train is

${y\text{?}\left( {t,1} \right)} = \left\{ \begin{matrix} {0.6t} & {0 \leq t \leq 5} \\ 3 & {5 < t \leq 20} \end{matrix} \right.$ ${y\text{?}\left( {t,2} \right)} = \left\{ \begin{matrix} {0.3t} & {0 \leq t \leq 10} \\ 3 & {10 < t \leq 20} \end{matrix} \right.$ ${y\text{?}\left( {t,3} \right)} = \left\{ \begin{matrix} {0.2t} & {0 \leq t \leq 15} \\ 3 & {15 < t \leq 20} \end{matrix} \right.$ ?indicates text missing or illegible when filed

and a reference trajectory signal is given by a waveform generator.

The formula (25) is solved to obtain the following PD-type ILC learning gain and observer gain of a nominal system:

$K_{1} = \begin{bmatrix} 0.403 & {- 40.08} & {- 1.952} & 0.003 & 0.036 & 0.002 \\ {- 1.951} & 0 & 3.944 & {- 39.969} & {- 1.993} & {- 0.002} \\ 0.036 & 0.002 & {- 1.993} & 0 & 1.549 & {- 40.019} \end{bmatrix}$ ${K_{2} = \begin{bmatrix} 24.578 & 0.029 & {- 0.002} \\ 0.029 & 25.746 & 0.001 \\ {- 0.002} & 0.003 & 24.692 \end{bmatrix}},{K_{3} = \begin{bmatrix} 24.538 & 0.03 & {- 0.002} \\ 0.03 & 25.706 & 0.002 \\ {- 0.002} & 0.003 & 24.652 \end{bmatrix}}$ $L = \begin{bmatrix} {- 0.009} & 0.024 & {- 0.006} & 0.002 & {- 0.006} & 0.007 \\ {- 0.006} & 0.002 & {- 0.01} & 0.029 & {- 0.006} & 0.002 \\ {- 0.006} & 0.007 & {- 0.006} & 0.002 & {- 0.009} & 0.024 \end{bmatrix}^{T}$

P-type ILC learning gain and corresponding observer gain are

$K_{1} = \begin{bmatrix} 0.504 & {- 40.1} & {- 1.995} & 0.008 & {- 0.002} & {- 0.001} \\ {- 1.995} & 0.012 & 1.788 & {- 40.119} & {- 2.006} & {- 0.012} \\ {- 0.002} & {- 0.001} & {- 2.006} & 0.007 & 0.371 & {- 40.108} \end{bmatrix}$ ${K\text{?}} = \begin{bmatrix} 25.284 & 0.006 & {- 0.026} \\ 0.006 & 25.238 & 0.006 \\ {- 0.026} & 0.006 & 25.268 \end{bmatrix}$ $L = \begin{bmatrix} {- 0.01} & 0.025 & {- 0.006} & 0.003 & {- 0.006} & 0.007 \\ {- 0.006} & 0.003 & {- 0.01} & 0.029 & {- 0.005} & 0.003 \\ {- 0.008} & 0.008 & {- 0.008} & 0.003 & {- 0.01} & 0.029 \end{bmatrix}^{T}$ ?indicates text missing or illegible when filed

The above iterative learning controller is implemented by an STM32F103RCT6 chip. An input signal of the chip is collected by a BENTLY 74712 speed sensor, stored into a stm32 chip through a conditioning circuit for calculation and used to form an iterative learning update law, and a signal obtained after calculation by a CPU (Central Processing Unit) is used as the control signal U_(k+1)(t) for a current trial. The control signal acts on a stepping motor DM3622 through a D/A (Digital/Analog) conversion circuit and is used for updating speeds of train particles until the given speed reference trajectory is reached.

FIG. 3 is an output curve of a first train particle under a nominal situation, FIG. 4 is an output curve of a second train particle under the nominal situation, and FIG. 5 is an output curve of a third train particle under the nominal situation. In order to further evaluate a tracking performance of the system, an RMS (Root-Mean-Square) error performance index is introduced:

${RMS} = {\sum\limits_{s = 1}^{3}\sqrt{\frac{1}{400}{\sum\limits_{i = 1}^{400}{e_{k}^{2}\left( {t,s} \right)}}}}$

FIG. 6 is an RMS contrast curve of a nominal system, and it is observed that with the presence of an input delay constant, a state of the system after T steps estimated by a state observer is used as feedback to advance output response of the system, so as to compensate input delay. With an increase of iterations, the control signal is updated constantly, and output of each train gradually reaches the desired speed trajectory, which proves effectiveness of the method according to the present invention. Furthermore, compared with P-type ILC, more tracking error information is used in PD-type ILC, and an RMS error converges faster along the trials, thus achieving a perfect tracking performance.

The formula (34) is solved to obtain PD-type ILC learning gain and observer gain of an uncertain system:

$K_{1} = \begin{bmatrix} 1.845 & {- 27.653} & {- 1.793} & {- 0.173} & {- 0.052} & {- 0.129} \\ {- 1.792} & {- 0.173} & 3.585 & {- 27.61} & {- 1.792} & {- 0.173} \\ {- 0.052} & {- 0.129} & {- 1.793} & {- 0.172} & 1.845 & {- 27.653} \end{bmatrix}$ ${K_{2} = \begin{bmatrix} 16.409 & 0.594 & {- 0.092} \\ 0.594 & 15.722 & 0.594 \\ {- 0.092} & 0.594 & 16.409 \end{bmatrix}},{K_{3} = \begin{bmatrix} 12.84 & 0.574 & {- 0.118} \\ 0.574 & 12.148 & 0.574 \\ {- 0.118} & 0.574 & 12.84 \end{bmatrix}}$ $L = \begin{bmatrix} {- 0.001} & 0.01 & 0.001 & {- 0.008} & {- 0.003} & 0.002 \\ 0.001 & {- 0.008} & {- 0.005} & 0.02 & 0.001 & {- 0.008} \\ {- 0.003} & 0.002 & 0.001 & {- 0.008} & {- 0.001} & 0.01 \end{bmatrix}^{T}$

Similarly, P-type ILC learning gain and corresponding observer gain are

$K_{1} = \begin{bmatrix} 1.851 & {- 27.954} & {- 1.799} & {- 0.231} & {- 0.052} & {- 0.139} \\ {- 1.799} & {- 0.231} & 3.598 & {- 27.861} & {- 1.799} & {- 0.231} \\ {- 0.052} & {- 0.138} & {- 1.799} & {- 0.231} & 1.852 & {- 27.953} \end{bmatrix}$ $K_{3} = \begin{bmatrix} 14.111 & 0.69 & 0.142 \\ 0.69 & 13.279 & 0.69 \\ {- 0.142} & 0.69 & 14.111 \end{bmatrix}$ $L = \begin{bmatrix} {- 0.001} & 0.008 & 0.001 & {- 0.008} & {- 0.002} & 0.002 \\ 0.001 & {- 0.008} & {- 0.004} & 0.018 & 0.001 & {- 0.008} \\ {- 0.002} & 0.002 & 0.001 & {- 0.008} & {- 0.001} & 0.008 \end{bmatrix}^{T}$

FIG. 7 is an output curve of the first train particle under an uncertain situation, FIG. 8 is an output curve of the second train particle under the uncertain situation, FIG. 9 is an output curve of the third train particle under the uncertain situation, and FIG. 10 is an RMS contrast curve of the uncertain system.

It is observed that with the simultaneous presence of the time-delay constant and structure uncertainty, a future state of the system is estimated by the state observer to form feedback acting on the system, such that the system outputs timely response after the time delay t, thus improving a control process. With the increase of the iterations, the output of each train asymptotically reaches the desired speed trajectory, and a tracking error converges along the trials, which shows the effectiveness and robustness for the structure uncertainty of the system of the method according to the present invention. In addition, almost 7 trials are required for perfect tracking by the PD-type ILC, and compared with the P-type ILC, convergence time is shorter, a convergence speed is higher, and the tracking performance is better.

The above description is only a preferred embodiment of the present application, and the present invention is not limited to the above embodiment. It is to be understood that other improvements and variations directly derived or suggested to those skilled in the art without departing from the spirit and idea of the present invention should be considered to be included within the protection scope of the present invention. 

What is claimed is:
 1. An ILC method for a multi-particle vehicle platoon driving system, comprising: step 1: establishing a spatially interconnected system model of the multi-particle vehicle platoon driving system; wherein a dynamic equation of the multi-particle vehicle platoon driving system is described as m{umlaut over (x)} _(i)({tilde over (t)})=k(x _(i+1)({tilde over (t)})−x _(i)({tilde over (t)}))−k(x _(i)({tilde over (t)})−x _(i−1)({tilde over (t)}))+u _(i)({tilde over (t)}−τ) y _(i)({tilde over (t)})=x _(i)({tilde over (t)})  (1)  wherein X_(i)({tilde over (t)}) represents a position of an ith train, u_(i)({tilde over (t)}) is control input and represents traction of the train, τ is a time-delay constant caused by signal transmission, y_(i)({tilde over (t)}) is measurement output and represents a speed of the train, m denotes mass of each train and has a variation range±ω _(m), and k is a spring coefficient and has a variation range±ω _(k); sampling time T is selected, and the equation (1) is approximately discretized using a finite difference method: that is, ${\left( {{\overset{.}{x}}_{i}\left( \overset{\sim}{t} \right)} \right)_{t,s} = \frac{{x_{1}\left( {{t + 1},s} \right)} - {x_{1}\left( {t,s} \right)}}{T}},{\left( {{\overset{¨}{x}}_{i}\left( \overset{\sim}{t} \right)} \right)_{t,s} = \frac{{x_{2}\left( {{t + 1},s} \right)} - {x_{2}\left( {t,s} \right)}}{T}}$  wherein t and s are discrete time and a train serial number respectively, and substitution of the above formula into the equation (1) gives a partial recurrence eqation: $\begin{matrix} {{x_{1}\left( {{t + 1},s} \right)} = {{x_{1}\left( {t,s} \right)} + {{Tx}_{2}\left( {t,s} \right)}}} & (2) \end{matrix}$ ${x_{2}\left( {{t + 1},s} \right)} = {{- 2\frac{kT}{m}{x_{1}\left( {t,s} \right)}} + {x_{2}\left( {t,s} \right)} + {\frac{kT}{m}{x_{1}\left( {t,{s + 1}} \right)}} + {\frac{kT}{m}{x_{1}\left( {t,{s - 1}} \right)}} + {\frac{T}{m}{u\left( {{t - \tau},s} \right)}}}$ y(t, s) = x₂(t, s) assuming that information transferred between the trains is their position information, that is, interconnected variables are set as: w₊(t,s)=w.(t,s)=x₁(t,s), v₊(t,s)=x₁(t,s−1), and v.(t,s)=x₁(t,s+1), and an output variable is set as: y(t,s)−x₂(t,s), the equation (2) is transformed into an uncertain spatially interconnected system model, i.e., $\begin{matrix} {\begin{bmatrix} {x\left( {{t + 1},s} \right)} \\ {w\left( {t,s} \right)} \\ {q\left( {t,s} \right)} \\ {y\left( {t,s} \right)} \end{bmatrix} = {\begin{bmatrix} A_{11}^{s} & A_{12}^{s} & B_{11}^{s} & B_{12}^{s} \\ A_{21}^{s} & 0 & 0 & 0 \\ C_{11}^{s} & C_{12}^{s} & D_{11}^{s} & D_{12}^{s} \\ C_{21}^{s} & 0 & 0 & 0 \end{bmatrix}\begin{bmatrix} {x\left( {t,s} \right)} \\ {v\left( {t,s} \right)} \\ {p\left( {t,s} \right)} \\ {u\left( {{t - \tau},s} \right)} \end{bmatrix}}} & (3) \end{matrix}$ wherein ${A_{11}^{s} = \begin{bmatrix} 1 & T \\ {- 2\frac{kT}{m}} & 1 \end{bmatrix}},{A_{12}^{s} = \begin{bmatrix} 0 & 0 \\ \frac{kT}{m} & \frac{kT}{m} \end{bmatrix}},{A_{21}^{s} = \begin{bmatrix} 1 & 0 \\ 1 & 0 \end{bmatrix}},{B_{11}^{s} = \begin{bmatrix} 0 & 0 \\ \frac{T}{m} & {- \frac{T}{m}} \end{bmatrix}},{B_{12}^{s} = \begin{bmatrix} 0 \\ \frac{T}{m} \end{bmatrix}}$ ${C_{11}^{s} = \begin{bmatrix} {- 2\varpi_{k}} & 0 \\ {- 2\frac{\varpi_{m}k}{m}} & 0 \end{bmatrix}},{C_{12}^{s} = \begin{bmatrix} \varpi_{k} & \varpi_{k} \\ \frac{\varpi_{m}k}{m} & \frac{\varpi_{m}k}{m} \end{bmatrix}},{C_{21}^{s} = \begin{matrix} \left\lbrack 0 \right. & {\left. 1 \right\rbrack,{D_{11}^{s} = \begin{bmatrix} 0 & 0 \\ \frac{\varpi_{m}}{m} & {- \frac{\varpi_{m}}{m}} \end{bmatrix}},{D_{12}^{s} = \begin{bmatrix} 0 \\ \frac{\varpi_{m}}{m} \end{bmatrix}}} \end{matrix}}$ in the above formula, x(t,s) ∈ R^(n), u(t,s) ∈ R^(q) and y(t,s) ∈ R^(m) represent state, input and output variables of an sth subsystem respectively; τ(0<τ<α) denotes the time-delay constant; w(t,s)=[w₊ ^(T)(t,s) w^(T)(t,s)]^(T) and v(t,s)=[v₊ ^(T)(t,S) v⁻ ^(T)(t,s)]^(T)represent interconnection between adjacent subsystems, and satisfy w ₊(s)=v ₊(s+1) w(s)=v(s−1)  (4)  boundary conditions thereof are set as: v₊(1)=w⁻(1)=0 and v⁻(n)=w₊(n)=0, and n is a number of the subsystems; p(t,s) denotes a pseudo-input channel of structure uncertainty, q(t,s) denotes a pseudo-output channel of the structure uncertainty, and p(t,s)θ_(s) q(t,s)  (5)  wherein an uncertainty block is: θ_(s)=diag{δ_(I), I_(r), . . . , δ_(I)I_(r)}, δ_(I) ∈ R, |δ_(i)|≤1, i=1, . . . ,f; δ_(i) describes a variation of a dynamic parameter of the system, and I_(r), is a unit matrix with a dimension of r_(f); step 2: transforming the spatially interconnected system model; wherein the model (3) is transformed into an equivalent one-dimensional dynamic model using a lifting technology, and lifting vectors are defined as follows: X(t)=[x(t,1)^(T), x(t,2)^(T), . . . , x(t,n)^(T)]^(T) V(t)=[v(t,1)^(T), v(t,2)^(T), . . . , v(t,n)^(T)]^(T) W(t)=[w(t,1)^(T), w(t,2)^(T), . . . , w(t,n)^(T)]^(T) P(t)=[p(t,1)^(T), p(t,2)^(T), . . . , p(t,n)^(T)]^(T) Q(t)=[q(t,1)^(T), q(t,2)^(T), . . . , q(t,n)^(T)]^(T) U(t−τ)=[u(t−τ,1)^(T), u(t−τ,2)^(T), . . . , u(t−τ,n)^(T)]^(T) Y(t)=[y(t,1)^(T), y(t,2)^(T), . . . , y(t,n)^(T)]^(T) then, the whole uncertain spatially interconnected system model can be equivalently described by the following model: $\begin{matrix} {\begin{bmatrix} {X\left( {t + 1} \right)} \\ {W(t)} \\ {Q(t)} \\ {Y(t)} \end{bmatrix} = {\begin{bmatrix} A_{11} & A_{12} & B_{11} & B_{12} \\ A_{21} & 0 & 0 & 0 \\ C_{11} & C_{12} & D_{11} & D_{12} \\ C_{21} & 0 & 0 & 0 \end{bmatrix}\begin{bmatrix} {X(t)} \\ {V(t)} \\ {P(t)} \\ {U\left( {t - \tau} \right)} \end{bmatrix}}} & (6) \end{matrix}$ wherein A₁₁ = diag{A₁₁¹, …, A₁₁^(n)}, A₁₂ = diag{A₁₂¹, …, A₁₂^(n)}, A₂₁ = diag{A₂₁¹, …, A₂₁^(n)} B₁₁ = diag{B₁₁¹, …, B₁₁^(n)}, B₁₂ = diag{B₁₂¹, …, B₁₂^(n)}, C₁₁ = diag{C₁₁¹, …, C₁₁^(n)} C₁₂ = diag{C₁₂¹, …, C₁₂^(n)}, C₂₁ = diag{C₂₁¹, …, C₂₁^(n)}, D₁₁ = diag{D₁₁¹, …, D₁₁^(n)} D₁₂ = diag{D₁₂¹, …, D₁₂^(n)} the model (6) contains the interconnected variables, and therefore, the model (6) is required to be further simplified; based on an interconnection characteristic of the formula (4) and the boundary conditions thereof, an equation relationship between the interconnected variables is obtained: W(t)=ηV(t)  (7)  wherein η is a permutation matrix unrelated to time t; the formula (7) is substituted into (6) to obtain V(t)=η⁻¹ A ₂₁ X(t)  (8)  and then, the formula (8) is substituted into (6) to eliminate the interconnected variables W(t) and V(t), so as to obtain the following equivalent uncertain model: $\begin{matrix} {\begin{bmatrix} {X\left( {t + 1} \right)} \\ {Q(t)} \\ {Y(t)} \end{bmatrix} = {\begin{bmatrix} {\overset{\_}{A}}_{11} & B_{11} & B_{12} \\ {\overset{\_}{C}}_{11} & D_{11} & D_{12} \\ C_{21} & 0 & 0 \end{bmatrix}\begin{bmatrix} {X(t)} \\ {P(t)} \\ {U\left( {t - \tau} \right)} \end{bmatrix}}} & (9) \end{matrix}$ wherein ${\overset{\_}{A}}_{11} = {A_{11} + {A_{12}\eta^{1}A_{21}}}$ ${\overset{\_}{C}}_{11} = {C_{11} + {C_{12}\eta^{- 1}A_{21}}}$ according to the formula (5), the following formula is obtained: P(t)=θQ(t)  (10)  wherein the uncertainty block is: θ=diag{θ₁, . . . ,θ_(n)}, θ_(i)≤I, i=1, . . . ,n; the formula (10) is substituted into (9) to obtain Q(t)=(I−D ₁₁θ)¹( C ₁₁ X(t)+D ₁₂ U(t−τ))  (11) then, the formula (11) is substituted into (9) to eliminate uncertainty variables P(t) and Q(t) using an elimination method, so as to obtain a general-form state space model: X(t+1)=(A+ΔA)X(t)+(B+ΔB)U(t−τ) Y(t)=CX(t)  (12)  wherein A=Ā₁₁, B=B₁₂, C=C₂₁, ΔA=B₁₁θ(1−D₁₁θ)⁻¹,C ₁₁ ΔB=B₁₁θ(I−D₁₁θ)⁻¹D₁₂ step 3: designing an ILC law based on a state observer; wherein the state space model (12) is described as an ILC structure form: X _(k+1)(t+1)=(A+ΔA)X _(k+1)(t)+(B+ΔB)U _(k+1)(t−τ) Y _(k+1)(t)=CX _(k+1)(t)  (13)  wherein k+1 represents a current operation trial of the system, t ∈ [0,α] denotes a limited working cycle of each trial of the system, and the input-delay constant satisfies a condition: τ<α; then, the ILC law is expressed as U _(k+1)(t)=U _(k)(t)+r _(k+1)(t)  (14)  the current control signal U_(k+1)(t) is equal to the control signal U_(k)(t) of a previous trial plus the update term r_(k+1)(t) which is calculated from previous error information; a system tracking error in a (k+1)th trial is e _(k−1)(t)=Y _(r)(t−τ)−Y _(k+1)(t)  (15)  wherein Y_(r)(t) is a desired output trajectory; considering delay in output response, the tracking error is re-described as e _(k−1)(t)=Y _(r)(t−τ)−Y _(k−1)(t)  (16) a state error vector is introduced: X _(k+1)(t+1)=X _(k+1)(t)−X _(k)(t)  (17) assuming that Y_(r)(0)=Y_(k)(0)=CX_(k)(0), and X _(k)(0)=0, that is, the system is returned to a same initial state in each trial, $\begin{matrix} {{{\overset{\_}{X}}_{k + 1}\left( {t + 1} \right)} = {{X_{k + 1}(t)} - {X_{k}(t)}}} & (18) \end{matrix}$  = (A + ΔA)[X_(k + 1)(t − 1) − X?(t − 1)] + (B + ΔB)[U_(k + 1)(t − 1 − τ) − U?(t − 1 − τ)] $= {{\left( {A + {\Delta A}} \right){{\overset{\_}{X}}_{k + 1}(t)}} + {\left( {B + {\Delta B}} \right){r_{k + 1}\left( {t - 1 - \tau} \right)}}}$ and $\begin{matrix} {{e_{k + 1}(t)} = {{Y_{r}\left( {t - \tau} \right)} - {Y_{k + 1}(t)}}} & (19) \end{matrix}$  = e_(k)(t) − C[X_(k + 1)(t) − X_(k)(t)] $= {{- {C\left( {A + {\Delta A}} \right)}{{\overset{\_}{X}}_{k + 1}(t)}} - {{C\left( {B + {\Delta B}} \right)}{r_{k + 1}\left( {t - 1 - \tau} \right)}} + {e_{k}(t)}}$ ?indicates text missing or illegible when filed in order to compensate input delay, the following state observer is constructed based on output information in the current trial: {tilde over (X)} _(k+1)(t+1)=A{tilde over (X)} _(k+1)(t)+Br _(k+1)(t−1)+L[ Y _(k+1)(t)− Y _(k+1)(t−τ)]  (20)  wherein Y _(k+1)(t)Y_(k+1)(t−1)−Y_(k)(t−1)CX _(k+1)(t), and {tilde over (X)}_(k+1)(t), is a τ-step ahead prediction of the state X _(k+1)(t); that is, {tilde over (X)}_(k+1)(t) is an estimation value of X _(k+1)(t+τ); and L is observer gain to be designed; an observation error is defined as {tilde over (e)} _(k+1)(t)= X _(k+1)(t)−{tilde over (X)} _(k+1)(t−τ)  (21) it is assumed that the update term in the ILC law (14) is r _(k+1)(r)=K ₁ {tilde over (X)} _(k+1)(t+1)+K ₂ e _(k)(t+τ)+K ₃ e _(k)(t+1+τ)−e _(k)(t+τ))  (22)  wherein K₁, K₂ and K₃ are learning gain to be designed; the update term is formed by state feedback information and PD-type previous tracking error information, and when the learning gain K₂ is equal to K₃, the formula (22) is simplified to P-type ILC; (22) is substituted into (20) to obtain $\begin{matrix} {{{\overset{\sim}{X}}_{k + 1}\left( {t + 1 - \tau} \right)} = {{\left( {A + {BK}_{1}} \right){{\overset{\sim}{X}}_{k + 1}\left( {t - \tau} \right)}} + {{LC}{{\overset{\sim}{e}}_{k + 1}\left( {t - \tau} \right)}} + {{B\left( {K_{2} - K_{3}} \right)}{e_{k}\left( {t - 1} \right)}} + {{BK}_{3}{e_{k}(t)}}}} & (23) \end{matrix}$ and $\begin{matrix} {{{\overset{\sim}{e}}_{k + 1}\left( {t + 1} \right)} = {{{\overset{\_}{X}}_{k + 1}\left( {t + 1} \right)} - {{\overset{\sim}{X}}_{k + 1}\left( {t + 1 - \tau} \right)}}} & (24) \end{matrix}$ $= {{\left( {{\Delta A} + {\Delta{BK}_{1}}} \right){{\overset{\sim}{X}}_{k + 1}\left( {t - \tau} \right)}} + {\left( {A + {\Delta A}} \right){{\overset{\sim}{e}}_{k + 1}(t)}} - {{LC}{{\overset{\sim}{e}}_{k + 1}\left( {t - \tau} \right)}} + {\Delta{B\left( {K_{2} - K_{3}} \right)}{e_{k}\left( {t - 1} \right)}} + {\Delta{BK}_{3}{e_{k}(t)}}}$ it is assumed that, ζ_(k+1)(t)=[X_(k+1) ¹(t−τ)e₁₊₁ ¹(t)]¹, and K=K₂-K₃, so as to obtain the following linear discrete repetitive process model with time delay: $\begin{matrix} {{\zeta_{k + 1}\left( {t + 1} \right)} = {{\hat{A}{\zeta_{k + 1}(t)}} + {{\hat{A}}_{t}{\zeta_{k + 1}\left( {t - \tau} \right)}} + {{\hat{B}}_{1}{e_{k}\left( {t - 1} \right)}} + {{\hat{B}}_{0}{e_{k}(t)}}}} & (25) \end{matrix}$ e_(k + 1)(t) = Ĉζ_(k + 1)(t) + D̂₁e_(k)(t − 1) + D̂₀e_(k)(t) wherein ${\hat{A} = \begin{bmatrix} {A + {BK}_{1}} & 0 \\ {{\Delta A} + {\Delta{BK}_{1}}} & {A + {\Delta A}} \end{bmatrix}},{{\hat{A}}_{t} = \begin{bmatrix} 0 & {LC} \\ 0 & {- {LC}} \end{bmatrix}}$ ${{\hat{B}}_{1} = \begin{bmatrix} {BK} \\ {\Delta{BK}} \end{bmatrix}},{{\hat{B}}_{0} = \begin{bmatrix} {BK}_{3} \\ {\Delta{BK}_{3}} \end{bmatrix}}$ Ĉ = [−C(A + ΔA) − C(B + ΔB)K₁ − C(A + ΔA)] D̂₁ = −C(B + ΔB)K, D̂₀ = I − C(B + ΔB)K₃ step 4: performing system stability analysis and learning gain solving on the linear discrete repetitive process model; wherein a Lyapunov function is selected as $\begin{matrix} {{V\left( {k,t} \right)} = {{V_{1}\left( {t,k} \right)} + {V_{2}\left( {k,t} \right)}}} & (26) \end{matrix}$ ${V_{1}\left( {t,k} \right)} = {{{\zeta_{k + 1}^{1}(t)}S{\zeta_{k + 1}(t)}} + {\sum\limits_{i = 1}^{t}{{\zeta_{k + 1}^{T}\left( {t - i} \right)}Q{\zeta_{k + 1}\left( {t - i} \right)}}}}$ V₂(k, t) = e_(k)^(T)(t − 1)P₁e_(k)(t − 1) + e_(k)^(T)(t)(P₂ − P₁)e_(k)(t) whereinS = diag{S₁, S₂} > 0, Q = diag{Q₁, Q₂} > 0, andP₂ > P₁ > 0;  V₁(t,k) represents an energy change along a trial, and V₂(k,t) represents an energy change between the trials; and increments of subfunctions are $\begin{matrix} {{\Delta{V_{1}\left( {t,k} \right)}} = {{{\zeta_{k + 1}^{T}\left( {t + 1} \right)}S{\zeta_{k + 1}\left( {t + 1} \right)}} - {{\zeta_{k + 1}^{T}(t)}\left( {S - Q} \right){\zeta_{k + 1}(t)}} - {{\zeta_{k + 1}^{T}\left( {t - \tau} \right)}Q{\zeta_{k + 1}\left( {t - \tau} \right)}}}} & (27) \end{matrix}$  = H_(k + 1)^(T)(t)Φ₁^(T)SΦ₁H_(k + 1)(t) − ζ_(k + 1)^(T)(t)(S − Q)ζ_(k + 1)(t) − ζ_(k + 1)^(T)(t − τ)Qζ_(k + 1)(t − τ) $\begin{matrix} {{\Delta{V_{2}\left( {k,t} \right)}} = {{{e_{k + 1}^{T}(t)}P_{2}{e_{k + 1}(t)}} - {{e_{k}^{T}\left( {t - 1} \right)}P_{1}{e_{k}\left( {t - 1} \right)}} - {{e_{k}^{T}(t)}\left( {P_{2} - P_{1}} \right){e_{k}(t)}}}} & (28) \end{matrix}$  = H_(k + 1)^(T)(t)Φ₂^(T)P₂Φ₂H_(k + 1)(t) − e_(k)^(T)(t − 1)P₁e_(k)(t − 1) − e_(k)^(T)(t)(P₂ − P₁)e_(k)(t) wherein ${{H_{k + 1}(t)} = \begin{bmatrix} {\zeta_{k + 1}(t)} \\ {\zeta_{k + 1}\left( {t - \tau} \right)} \\ {e_{k}\left( {t - 1} \right)} \\ {e_{k}(t)} \end{bmatrix}},{\Phi_{1} = \left\lbrack \begin{matrix} \hat{A} & {\hat{A}}_{\tau} & {\hat{B}}_{1} & {\left. {\hat{B}}_{0} \right\rbrack,{\Phi_{2} =}} & \left\lbrack \hat{C} \right. & 0 & {\hat{D}}_{1} & \left. {\hat{D}}_{0} \right\rbrack \end{matrix} \right.}$ a total increment of the function is $\begin{matrix} {{\Delta{V\left( {k,t} \right)}} = {{{\Delta{V_{1}\left( {t,k} \right)}} + {\Delta{V_{2}\left( {k,t} \right)}}} = {{H_{k + 1}^{T}(t)}{\prod{H_{k + 1}(t)}}}}} & (29) \end{matrix}$ wherein $\prod{= {{{\begin{bmatrix} \hat{A} & {\hat{A}}_{\tau} & {\hat{B}}_{1} & {\hat{B}}_{0} \\ \hat{C} & 0 & {\hat{D}}_{1} & {\hat{D}}_{0} \end{bmatrix}^{T}\begin{bmatrix} S & 0 \\ * & P_{2} \end{bmatrix}}\begin{bmatrix} \hat{A} & {\hat{A}}_{\tau} & {\hat{B}}_{1} & {\hat{B}}_{0} \\ \hat{C} & 0 & {\hat{D}}_{1} & {\hat{D}}_{0} \end{bmatrix}} - \begin{bmatrix} {S - Q} & 0 & 0 & 0 \\ * & Q & 0 & 0 \\ * & * & P_{1} & 0 \\ * & * & * & {P_{2} - P_{1}} \end{bmatrix}}}$ if ΔV(k,t)<0 is satisfied for all k and t, the model (25) is stable along the trials, which has an equivalent condition of Π<0; a Schur complement lemma is applied to Π<0, left and right sides of an inequality are multiplied by diag{I,I,S¹,S¹,P₂ ¹,P₂ ¹}, variable substitution is performed, and it is assumed that S₁ ⁻¹=W₁, S₂ ¹=W₂, S₁ ¹Q₁S₁ ⁻¹=X₁, S₂ ¹Q₂S₂=X₂, P₂ ¹P₁P₂ ¹=Z₁, P₂ ¹=Z₂, K₁S₁ ¹=K₁W₁=R₁, LCS₂ ¹=LCW₂=R, KP₂ ¹=K₂=R₂, K₁P₂ ¹=K₃Z₂=R₁, so as to obtain the following conclusion: for the nominal linear discrete repetitive process model with time delay described in the formula (25), if there exist matrices W=diag{W₁,W₂}>0, X=diag{X₁,X₂}>0, Z₁>0, Z₂>0 and R, R₁, R₂ and R₃ satisfying the following linear matrix inequality: $\begin{matrix} {\Xi_{1} = \text{ }{\begin{bmatrix} {- W_{1}} & 0 & 0 & {{AW}_{1} + {BR}_{1}} & 0 & 0 & R & {{BR}\text{?}} & {BR}_{1} \\ * & {- W_{2}} & 0 & 0 & {AW}_{2} & 0 & {- R} & 0 & 0 \\ * & * & {- Z_{2}} & {- {C\left( {{AW}_{1} + {BR}_{1}} \right)}} & {- {CAW}_{2}} & 0 & 0 & {- {CBR}_{2}} & {Z_{2} - {{CBR}\text{?}}} \\ * & * & * & {{- W_{1}} + X_{1}} & 0 & 0 & 0 & 0 & 0 \\ * & * & * & * & {{- W_{2}} + X_{2}} & 0 & 0 & 0 & 0 \\ * & * & * & * & * & {- X_{1}} & 0 & 0 & 0 \\ * & * & * & * & * & * & {- X_{2}} & 0 & 0 \\ * & * & * & * & * & * & * & {- Z_{1}} & 0 \\ * & * & * & * & * & * & * & 0 & {{- Z_{2}} + Z_{1}} \end{bmatrix} < 0}} & (30) \end{matrix}$ ?indicates text missing or illegible when filed  the model (25) is stable along the trials, and the learning gain of the update term (22) and the gain of the state observer (20) are: Ξ₁ +MΘN+N ^(T)Θ^(T) M ^(T)<0  (32)  when the structure uncertainty of the system is considered, the Schur complement lemma is applied to Π<0, the left and right sides of the inequality are multiplied by diag|I, I, S¹, S¹, P₂ ¹, P₂ ¹|,variable substitution is performed, and it is assumed that S₁ ¹=W₁, S₂ ¹=W₂, S₁ ¹Q₁S₁ ¹=X₁, S₂ ¹Q₂S₂ ¹=X₂, P₂ ¹P₁P₂ ¹=Z₁, P₂ ¹=Z₂, K₁S₁=K₁W₁=R₁, LCS₂ ¹=LCW₂=R, KP₂ ¹=KZ₂=R₂, K₃P₂ ¹=K₃Z₂=R₃, so as to obtain Ξ₁ +MΘN+N ^(T)Θ^(T) M ^(T)<0  (32)  wherein M=[0 B₁₁ ^(T)−B₁₁ ^(T)C¹ 0 0 0 0 0 0]^(r) N=[0 0 0 C ₁₁W₁+D₁₂R₁ C ₁₁W₂ 0 0 D₁₂R₂ D₁₂R₁] Θ=θ(I−D₁₁θ)⁻¹, θ^(T)θ≤I the formula (32) is equivalent to $\begin{matrix} {{\Xi_{1} + {{\begin{bmatrix} \varepsilon^{- 1} & N^{T} & {\varepsilon M} \end{bmatrix}\begin{bmatrix} I & {- D_{11}} \\ * & I \end{bmatrix}}^{1}\begin{bmatrix} {\varepsilon^{- 1}N} \\ {\varepsilon M^{T}} \end{bmatrix}}} < 0} & (33) \end{matrix}$ whereinε > 0; according to the Schur complement lemma, the formula (33) is described as $\begin{bmatrix} \Xi_{1} & {\varepsilon^{1}N^{T}} & {\varepsilon M} \\ * & {- I} & D_{11} \\ * & * & {- I} \end{bmatrix} < 0$ left and right sides of the above formula are multiplied by diag{I, εI, εI} and ε² is replaced by ε, so as to obtain the following conclusion: for the uncertain linear discrete repetitive process model with time delay described in the formula (25), if there exist the matrices W=diag{W₁,W₂}>0, X=diag{X₁,X₂}>0, Z₁>0, Z₂>0 and R, R₁, R₂ and R₃ satisfying the following linear matrix inequality: $\begin{bmatrix} \Xi_{1} & \Xi_{2} \\ * & \Xi_{3} \end{bmatrix} < 0$ wherein $\Xi_{2} = \begin{bmatrix} 0 & 0 & 0 & {{{\overset{\_}{C}}_{11}W_{1}} + {D_{12}R_{1}}} & {{\overset{\_}{C}}_{11}W_{2}} & 0 & 0 & {D_{12}R_{2}} & {D_{12}R_{3}} \\ 0 & {\varepsilon B_{11}^{T}} & {- \varepsilon B_{11}^{T}C^{T}} & 0 & 0 & 0 & 0 & 0 & 0 \end{bmatrix}^{T}$ $\Xi_{3} = \begin{bmatrix} {- \varepsilon I} & {\varepsilon D_{11}} \\ * & {- \varepsilon I} \end{bmatrix}$  the model (25) is robustly stable along the trials, and the learning gain of the update term (22) and the gain of the state observer (20) are given by the formula (31); step 5: giving specific parameters of the multi-particle vehicle platoon driving system, determining the learning gain of the ILC law and the corresponding observer gain, applying the control signal of the current trial to the ILC state space model to obtain the output of the current trial, and then repeatedly adjusting and controlling the ILC law to allow the output of each train of the vehicle platoon driving system to asymptotically reach the desired speed trajectory. 